Predictors of drought-induced crop yield/losses in two agroecologies of southern Tigray, Northern Ethiopia

The consequences of prolonged precipitation-deficient periods are primarily substantial water deficit. The spatial characteristics of drylands and various socioeconomic factors worsen droughts’ impacts and deepen poverty among agrarian communities, with attendant food security (stability dimension) implications. This study utilizes a combination of climate, remote sensing and field survey data to obtain first-hand information on the impacts of recent (2015 and 2017) droughts on crop yield in southern Tigray, northern Ethiopia. Annual and seasonal rainfall, annual and seasonal Normalized Difference Vegetation Index (NDVI) and Deviation of NDVI (Dev-NDVI), and monthly Standardized Precipitation Index (SPI) (SPI-1, SPI-3 and SPI-12) for June to October, were considered as likely factors that could relate with yield and yield loss in the area. Correlation and multiple linear stepwise regression statistical techniques were used to determine drought-yield relationships, and identify more accurate predictors of yield and yield losses in each of the drought years. The area witnessed a more widespread precipitation deficit in 2015 than in 2017, where the lowland area recorded entire crop (sorghum) losses. Also, droughts manifested spatiotemporal variations and impacts across the two different agroecologies—primarily reduction in vegetation amounts, coinciding with the planting and maturing stages of barley and sorghum. Crop failures, therefore, translated to food shortages and reduced income of smallholder farmers, which denotes food insecurity in the time of droughts. Seasonal rainfall and June Dev-NDVI predicted 66.9% of 2015 barley and sorghum yield-loss, while NDVI predicted 2017 sorghum yield by 96%. Spate irrigation should be further popularized in the low-lying areas of Raya Azebo to augment for future deficiencies in the kiremt rainfall.

Global climate change has resulted in the increasing occurrence of extreme weather and climatic events such as floods and droughts. According to Eze et al. 1 and Negash et al. 2 , outputs of rain-fed agriculture are at risk of influence by drought events. Droughts occur when a shortage of precipitation is received in a given season, or when there is a departure from a known average. An abnormally prolonged dry condition is designated as drought by Bayissa et al. 3 , whereas in the view of Wilhite 4 , the timing and the effectiveness of precipitation could be linked to drought, hence making each drought event unique in impacts and characteristics. The consequences of prolonged precipitation-deficient periods are primarily substantial water deficit. Water shortages are exacerbated by the spatial characteristics of drylands, especially in arid and semi-arid environments. Additional disadvantages such as a largely rural population with very little education, totally dependent on agriculture for living and livelihood, and a limited capacity to withstand droughts could further extend impacts of droughts.

Materials and methods
The study area. This study was implemented in two purposively selected woredas (districts) in the southern Tigray region, northern Ethiopia. Tigray has been previously established as one of the most drought-vulnerable regions in Ethiopia as described in earlier sections of this paper. The study area is located between 12° 38′ 44″-12° 57′ 10″ N Latitude and 39° 27′ 18″-39° 55′ 56″ E Longitude, with an altitude ranging from 1109 to 3760 m above sea level (m.a.s.l.). The two purposively-selected woredas (districts) were to represent the highland and lowland agroecology for this study. The study area thus includes one district of mid and highland agroecology (Endamehoni) and another woreda (Raya Azebo) of lowland agroecology. Nine (9) tabias/kabeles (villages) are included in the portion selected for this study (Fig. 1).

Data sources and analyses methods. Meteorological data and meteorological drought indicator. The
Climate Hazards group InfraRed Precipitation with Stations (CHIRPS) grid precipitation dataset from 2015 to 2018, obtained from (ftp:// ftp. chg. ucsb. edu/ pub/ org/ chg/ produ cts/ CHIRPS-2.0/) was used for this study. The CHIRPS is the combination of a "high-resolution climatology, time-varying cold cloud duration precipitation estimates, and in situ precipitation estimates", with high research quality and very low errors 19,20 . The high reliability, spatial spread and accuracy of the CHIRPS over the point-based data makes it preferred for this study. Also, The CHIRPS has been described as dependable and accurate for use in Ethiopia, especially for obtaining data for areas without meteorological stations 21,22 . Furthermore, despite the resolution of the CHIRPS dataset being at 5 km spatial resolution, Dinku et al. 23 conducted a validation of CHIRPS across eastern Africa and reports that the dataset offers a higher skill than other alternative rainfall products with low or no bias for Ethiopia. Moreover, Peng et al. 24  www.nature.com/scientificreports/ studies, hence a justification for use in this study. The rainfall variables used in this study include the annual and main (kiremt) season (JJAS) rainfall. This data is used to assess meteorological drought using the Standardized Precipitation Index (SPI). The commonly used meteorological drought indicators include the Standardized Precipitation Index (SPI) and the Standardized Precipitation Evapotranspiration Index (SPEI). The SPI is recommended by the World Meteorological Organization 25 to be used by all National Meteorological and Hydrological Services around the world to characterize meteorological droughts. Designed by McKee et al. 26 to quantify the magnitude of dry and wet conditions in the desired location, the SPI is effective in drought detection in multiple time-scales, generally between 1 and 48 months. The WMO and GWP 27 posit that SPI values for 3 months (SPI-3) or less, might be useful for basic drought monitoring; while values for 6 months (SPI-6) or less, are fit for monitoring agricultural impacts; and values for 12 months (SPI-12) or longer, are ideal for hydrological impacts. Besides the abovementioned points, Yildirak and Selcuk-Kestel 28 regard SPI as easy to implement and ideal for the determination of drought risks due to its 'intrinsic probabilistic nature' and favourability in making spatiotemporal comparisons among diverse climatic condition. These inform the preference for SPI in this study for meteorological drought assessment. We however reiterate the limitation that since cumulative precipitation is the sole input, SPI is deficient when accounting for the temperature component of droughts, as water balance/use are not accounted for.   29 , MODIS data products are accurate and appropriate for drought monitoring. Following download, the images were re-projected and resampled from a sinusoidal to geographic projection (WGS84); and then NDVI values were extracted from the re-projected (*.tif) images using R scripts. To obtain the conventional NDVI values ranging from − 1 to + 1, the extracted values were rescaled with the factor 0.0001. The NDVI variables used in this study include the Annual NDVI, the Dev-NDVI of the main cropping (kiremt) season, and the monthly Dev-NDVI of June to October. A lot of satellite-based indicators have been developed, with the ability to accurately detect droughts over large spatiotemporal resolutions. The Normalized Difference Vegetation Index (NDVI) is a popular droughtmonitoring index developed by Kogan 30 , for determination of droughts by the measurement of vegetation vigor and cover in an index ranging from − 1 to + 1. However, the NDVI itself needs to be computed vis-à-vis anomalies to indicate the presence of drought and its severity. The Dev-NDVI, which is the deviation of NDVI from its long-term mean will easily indicate the magnitude of wetness or dryness of an area. Hence, a negative Dev-NDVI is an indicator of below-normal vegetation condition/health, thus suggesting a drought situation. The magnitude of a drought spell is defined by the degree of negative departure from the long-term mean NDVI 31,32 . Kourouma et al. 33 assert that 'vegetation in general, and crops productivity, in particular, are responsive to change in NDVI values' . The NDVI and Dev-NDVI are therefore used in this study to obtain a clear perception of agricultural droughts.
Farmers' yield (survey) data. Self-reported farmers' production data were obtained from a survey of 34 smallholders for this study. A minimum N ≥ 25 is required for more variance and is recommended for regressions or meta-regressions 34 . Hence, a sample size of 34 is considered optimal for our study. A one-time direct interview of 34 farmers was conducted to elicit information on the dry grain weight of crop yield measured in quintals per hectare. After maturity, farmers harvest their grains and allow them dry in a part of the farmland. The grains are beaten for separation from their husks and then bagged. Farmers are able to weigh their grains before storage and sales.
Farm owners used in this study (sorghum farmers = 16; barley farmers = 18) were randomly selected, and are those with high yield recall from memory. Ethical procedures such as obtaining their consent to participate in the study and assurance of anonymity were strictly followed for all participants in the study. While more farmers were interviewed as they are always found in groups, only those who met our criteria of crop type and recent crop yield memory were included in the study. Their farmlands were measured with the use of HP-IPAQ (a pocket PC enabled with GIS features for plotting). The device was useful for appropriate relation of the yield to specific farmlands. Two cereal crops-barley and sorghum were selected in the highland and lowland respectively, due to their widespread importance as dominant crops and staple foods in the area, and also to improve the accuracy of yield recall from farmers' memory.
The village (Tabia) or plot-level data are not available at either the Agricultural Bureau or the Central Statistics Agency. The use of district-level data available in these offices would not be representative of the local situations. This explains the choice of plot-level data in the study area, which is considered of optimum quality. Eze et al. 1 adopted the same data to successfully develop proportionally-derived area-yield insurance index with high success in estimating payouts even for those who reported 100% losses. Also, Eze et al. 35 attempted to develop remote crop yield estimation using crop water requirements using this data. From the evaluation of the models, there is a good fit with the survey data with a coefficient of determination (R 2 ) of up to 88%. This data is thus applied in this study with the confidence of optimal performance. All methods in this study were carried out in accordance with relevant guidelines and regulations for conducting surveys. In particular, informed consent was obtained from all human subjects, and anonymity has been maintained in presenting the results obtained from analyses of collected data. We confirm that the academic board of the Institute of Climate and Society, Mekelle University, Ethiopia orally approved the conduct of this study during the presentation of the proposed research, survey instrument and method for data collection. Neither plants nor animals were used in this study, which ruled out the need for written ethical clearance or approval. Data analyses. Data pre-processing procedures. The analysis of historical meteorological droughts was done using the Standardized Precipitation Index (SPI) in R software version 3.4 with the aid of suitable scripts and packages. The SPI was calculated on three timescales: 1-, 3-and 12-month (SPI-1, SPI-3, and SPI-12). To obtain agricultural droughts Dev-NDVI was calculated using the Eq. (1).
where the NDVI i is the NDVI value for month 'i'; NDVI mean-m is the long-term mean NDVI for the same month m over the spatiotemporal period.
Three operations were implemented to get three variables of yield data for the analyses.
1. Using the raw yield data obtained from farmers, conversion to yield per hectare was done by dividing the actual (raw) yield value by the area of land and multiplying the result by 10,000 (which is the area of a hectare). Hence the formula: where Yj is the crop yield in j Year of a Tabia (village); Ȳ is the average, and σ is the standard deviation of crop yield during the period under review.

Correlating yield, climatic and drought-indicator variables.
To obtain the relationship between drought and crop yield in the study area, YpHa, SCY and YLR for 2015 and 2017 were correlated with annual rainfall; kiremt rainfall; annual NDVI; annual Dev-NDVI; kiremt Dev-NDVI; monthly Dev-NDVI (for June to October); and monthly SPI-1, SPI-3 and SPI-12 (for June to October). The use of all these variables in measuring the relationship is to determine which factor is most related to the declared farmers' yield and losses for the years under review. The analysis of the relationship between variables listed above with YpHa, SCY and YLR was done in three batches. Firstly, the relationships across the entire study area were assessed. Next, a split Woreda-based analysis was conducted to ascertain the relationships based on the two different woredas (Endamehoni and Raya Azebo) selected for the study. The purpose of the split-site analysis is to account for inherent variations expected in rainfall patterns, crop yield and crop loss in the two agroecologically dissimilar districts (woredas).
The non-parametric (Spearman's) correlation analysis was implemented, as this would not demand bivariate normality of the dataset. Some previous studies 37-40 all adopted the Spearman's correlation for their analyses. Interpretation of the strength of association among variables as reflected by the resultant correlation coefficients adopted in this study was obtained from the groupings of Cohen et al. 41 as shown below (Table 1).

Regressing yield, climatic and drought-indicator variables.
Multiple linear stepwise regression was implemented in SPSS to obtain the most significant and relevant predictors for both crop yield and yield loss in the drought years of 2015 and 2017. The stepwise regression is an approach to selecting a subset of effects for a regression model. This approach has been widely used in literature, with several authors preferring its application where little theory to guide the selection of terms for a model exists; or where there is the focus on predictors which interactively seem to provide a good fit; and where the performance of a model is to be improved by reducing the variance from unnecessary terms [42][43][44][45][46] . The stepwise regression also exempts variables of high multicollinearity (variance inflation factor-VIF) in the selection of the most relevant predictors, hence its adoption in this study. The step-wise multiple regression was therefore conducted using the threshold of probability of the predictor to enter as ≤ 0.050 and the probability of the predictor to be removed as ≥ 0.100 (Table 2).

drought and crop yield: relationship and climatic/vegetative predictors. The results
are presented in three subsections to represent the lowland case study (Raya Azebo), the highland case study (Endamehoni), and the entire study area. These sections correspond to the three-batch analyses implemented to unravel the relationship of drought with crop yield and arriving at key drought indicators that more accurately predicted the crop yield in the event of drought. Figure 2 contains the key spatial variables used to obtain the results.
Area m 2 × 10, 000  (Table 3). This implies that higher annual rainfall translates to higher crop yield while higher meteorological drought in June and September led to lower crop yield. Kiremt season rainfall also had a strong positive significant relationship with YpHa. Other notable variables indicating a moderate positive relationship with YpHa and SCY are annual Dev-NDVI, kiremt Dev-NDVI; while a moderate negative relationship exists between July and August SPI-1, June and July SPI-1; June and October SPI-3 (Table 3). Annual rainfall and kiremt season rainfall had a strong negative relationship with YLR. Hence, the lower the (annual and kiremt)  www.nature.com/scientificreports/ rainfall, the higher the yield loss. On the other hand, annual Dev-NDVI, June and October SPI-1 are strongly positively and significantly associated with YLR. Hence, a higher deviation from NDVI and higher meteorological drought reflected in higher crop losses. Kiremt Dev-NDVI, June Dev-NDVI, July SPI-1, June and October SPI-3, and June SPI-12 were moderately positively correlated with YLR. Annual rainfall and July SPI-3 are the most significant predictors of SCY in Raya Azebo. The two-step regression (Supplementary Table S.01) results show that annual rainfall alone explains 38% of the variation in the standardized crop yield while and the combination of both predictors-annual rainfall and July SPI-3 for the year explain 56.3% of the variation in SCY for Raya Azebo (Table 5). Similar to the SCY, the YpHa is significantly predicted by annual rainfall and July SPI-3 with the coefficient of determination (r 2 ) indicating that 35.7% of the variations in the yield per hectare obtained in the study area is attributed to annual rainfall alone, while the combination of annual rainfall and July SPI-3 explains 53.3% of the variations (Table 5). For the 2015 season, about 75% of farmers sampled in Raya Azebo had zero yields and 100% of crop losses. The stepwise regression for YLR was unsuccessful due to the multiple zero figures.
Highland (Endamehoni). There are slightly different results obtained in the analyses of data in this area. There is a strong negative significant relationship between kiremt Dev-NDVI, June Dev-NDVI, July and September SPI-1 and SPI-3; and July September and October SPI-12, all with YpHa ( Table 4). The SCY had a strong negative significant relationship with kiremt Dev-NDVI, September SPI-1; July and September SPI-3; and July, September and October SPI-12. Additionally, it is worthy of mention that the correlation between kiremt season rainfall, August SPI-1, June, August and October SPI-3, and August SPI-12 with YpHa and SCY was moderately negative, though not significant. Therefore, higher wetness led to lower yield in this woreda. The crop yield loss was strongly and positively significantly related with Annual Dev-NDVI, June Dev-NDVI, August and September SPI-1; June to October SPI-3, and August to October SPI-12. Furthermore, kiremt season rainfall, July SPI-1 and July SPI-12 showed a moderate positive relationship with the YLR but not significant. Higher rainfall reflected in higher yield loss in this area and higher Dev-NDVI did not imply higher yield as crop losses were recorded due to excessive rainfall, as confirmed by data obtained during farmers' survey (Table 4).  Table S.01) indicates September SPI-12 and October Dev-NDVI as the key predictors of yield loss. The coefficient of determination (r 2 ) show that September SPI-12 alone explains 38.4% of the variation in yield loss ratio; while the combination of September SPI-12 and October Dev-NDVI explain 62.4% (Table 5).
Entire study area. From the results obtained from the entire area analysis, there is a strong positive significant relationship between Annual Rainfall and Kiremt Rainfall with farmers' YpHa and SCY (Table 2). This implies that higher rainfall led to higher crop yield. Conversely, there is a strong negative significant relationship between annual rainfall, kiremt rainfall and farmers' YLR. Besides, the annual NDVI shows a strong positive significant relationship between YpHa and SCY while annual Dev-NDVI shows a moderate/medium positive significant relationship between the same variables. This implies that the deviation of vegetation index from the long-term average does not always translate to yield reduction. Farmers' YLR shows a moderate negative significant relationship with Annual NDVI and a weak negative relationship with Annual Dev-NDVI (Table 2).
Furthermore, the monthly Dev-NDVI indicate various relationships with farmers' YpHa, SCY and YLR. The SPI at the 1, 3, and 12-month scales show different levels of relationship for different months with the YpHa, SCY and YLR. At the 1-month scale of SPI, there is a strong negative significant relationship between YpHa, SCY and SPI-1 values for June and October, with August showing a nearly strong negative significant association (Table 2). This implies that higher SPI-1 drought during the planting and maturing stages impacted crops' performance and resulted in yield reduction. A strong positive significant relationship exists between YLR and SPI-1 for June, August and October meaning that higher droughts in these months, which are crucial growing and maturing phases of the crops studied, translates to higher yield losses. The 3-month scale of SPI shows a strong positive significant relationship with YpHa and SCY in September and a strong negative significant relationship with YLR. For the 12-month scale of SPI, a strong positive significant relationship exists between YpHa and SCY in July, September and October (Table 2). While the YLR has a strong negative significant relationship with July and September SPI-12, it has a medium negative significant relationship with October SPI-12 values (Table 2). Hence, higher SPI-12 drought in July and September led to higher crop losses, whereas October drought had no major impact on crop losses.
Kiremt season rainfall, June Dev-NDVI and annual NDVI are the key predictors for SCY in the entire study area. From the three-step regression, kiremt rainfall alone explains 39.8% of the variation in the standardized crop yield; Kiremt rainfall and June Dev-NDVI explain 50.4%; and the combination of all three predictors-kiremt rainfall, June Dev-NDVI and the annual NDVI for the year explain 57.5% of the variation in crop yield in the study area. The key predictors for YLR are kiremt rainfall and June Dev-NDVI. The two-step regression shows that kiremt rainfall alone explains 51.3% of the variation in yield loss ratio; while the combination of kiremt rainfall and June Dev-NDVI explain 66.9% (Table 5). Some additional variables that are excluded from the step-wise threshold, though with significant predictive possibilities but high multicollinearity, include the September SPI-3 (p = 0.047) and October SPI-12 (p = 0.046). From a one-step regression analysis (Supplementary Table S.01), annual rainfall is the most significant predictor of YpHa. The coefficient of determination (r 2 ) obtained shows that 40.3% of the variations in the yield per hectare obtained in the entire study area is attributed to annual Rainfall. Table 5. Predictors of YpHa, SCY and YLR in the 2015 drought year in the study area. ***Regression predictor is significant at the above 0.001 level. **Regression predictor is significant at the 0.01 level. *Regression predictor is significant at the 0.05 level.  (Fig. 3) is inconsistent with the 2015 drought year (Fig. 2) as no location recorded 100% crop loss. Also, areas that had higher crop losses in the previous drought of 2015 had better and improved yield in this season. As earlier indicated above, the results are presented in three subsections to represent the lowland case study (Raya Azebo), the highland case study (Endamehoni), and the entire study area. These sections also correspond to the three-batch analyses implemented to unravel the relationship of drought with crop yield and arriving at key drought indicators that more accurately predicted the crop yield in the drought events (Table 6).

Raya Azebo (lowland).
There is a strong positive significant relationship between annual NDVI with YpHa and SCY. A positive moderate not-significant relationship also exists between kiremt Dev-NDVI, YpHa and SCY (Table 7). June to September Dev-NDVI values show a moderate negative not-significant relationship with YpHa  www.nature.com/scientificreports/ and SCY. Crop loss relates moderately with kiremt season rainfall, June and July Dev-NDVI, August SPI-1; July to September SPI-12 values. All these were not significant in their relationship. However, June SPI-3 has a strong negative significant relationship with YLR while October SPI-12 has a strong positive relationship with YLR, which is also significant ( Table 7). The results of a seven-step multiple regression (Supplementary Table S.02) indicate that annual NDVI for 2017, and June to October Dev-NDVI values are key predictors of SCY. The combination of annual NDVI, June Dev-NDVI, August Dev-NDVI, July Dev-NDVI and October Dev-NDVI explains 96% of the variation of SCY in Raya Azebo for this farming season. Similarly, a seven-step regression presents the annual NDVI and Dev-NDVI, June SPI-12, June, August and September Dev-NDVI as the most significant predictors of YpHa with 96.6% explanation of the variation in yield per hectare obtained in Raya Azebo ( Table 9). The key predictor of crop losses in Raya Azebo during the 2017 drought was Kiremt Rainfall ( Table 9). The coefficient of determination (r 2 ) obtained shows that 37.1% of the variations in the yield loss ratio obtained in the study area in 2017 is attributed to the kiremt rainfall values (Table 9).
Endamehoni (highland). Key results found in this area indicate that YpHa and SCY are moderately positively correlated with Annual and June Dev-NDVI (Table 8). The SCY, however, shows a strong positive relationship with kiremt Dev-NDVI, August and September Dev-NDVI. The YLR had a strong negative significant relationship with kiremt Dev-NDVI, September SPI-3; July, September and October SPI-12. The variable also had a moderate negative correlation with annual and June Dev-NDVI; July, September and October SPI-1; October SPI-3 values.
The most relevant predictors of SCY in Endamehoni are September Dev-NDVI and the 2017 Kiremt Dev-NDVI in a two-step computation (Supplementary Table S.02). The coefficient of determination (r 2 ) shows that only the September Dev-NDVI is a predictor of SCY explain 29.5% of the variation, while the addition of the kiremt season Dev-NDVI for the year explains 46.9% variation in standardized crop yield in the area for the year under review (Table 9). Similarly, the September Dev-NDVI is the sole predictor of YpHa in the area with www.nature.com/scientificreports/ 45.2% of variations explained in the one-step regression implemented (Table 9). September SPI-3 is the only, most relevant and significant predictor of crop losses in this location for the year 2017. The coefficient of determination (r 2 ) shows that September SPI-3 explains 45.4% of the variation in yield loss ratio in the area (Table 9). A nearly-significant predictor is the June Dev-NDVI with (p = 0.052).
Entire study area. The results depict a moderate positive significant relationship between annual rainfall, kiremt season rainfall, annual NDVI and YpHa (Table 6). Similarly, a moderate positive significant relationship exists between Annual NDVI and SCY. Also, the July, September and October SPI-1; September and October SPI-3; together with June, July and October SPI-12 all show a moderate positive relationship with YpHa and SCY. In correlating crop loss with the various factors, only the annual rainfall and kiremt season rainfall showed a positive significant relationship with YLR, while June SPI-3 values show a moderate positive significant relationship. Annual NDVI and July SPI-12 are the key predictors of SCY from the two-step multiple regression analysis implemented (Supplementary Table S.02). Annual NDVI alone explains 19.6% of the variation in the standardized crop yield while the combination of both predictors-annual NDVI and the July SPI-12 values for the year explain 33.1% of the variation in SCY (Table 9). However, for the YpHa, July SPI-12 as the only and most relevant and significant predictor with the coefficient of determination (r 2 ) obtained shows that 24.4% of the variations in the yield per hectare obtained in the study area is attributed to the July SPI-12 values ( Table 9).
The most significant predictors of yield loss in the study area are kiremt rainfall and September SPI-3. The coefficient of determination (r 2 ) show that kiremt rainfall alone explains 19% of the variation in yield loss ratio; while the combination of kiremt rainfall and September SPI-3 explain 33% variations in yield loss within the entire study area (Table 9).

Discussion
The year 2015 is adjudged as the most recent drought year of high magnitude and one of the worst events over the past 30 years covering a large area and exacerbating food insecurity in Ethiopia 3,47 . These literature reports are supported by evidence of low crop yield in a larger part of the study area (Fig. 2). Also, the year 2017 is recorded as the most recent drought year in some parts of Tigray/Ethiopia with Gross Domestic Product (GDP) losses due to climate variability estimated to be around 1-4% excluding human losses 8 . These events were expected to have impacted on crop yield and yield losses in the study area.
From this study, drought has again been proven to be location-specific, and not a broad spatio-temporal phenomenon at any given time. In an earlier study, Viste et al. 5 describe most historic droughts as being more of 'local or regional character' with dissimilar effects at different places and seasons. Each drought year reviewed (2015 and 2017) stands out uniquely in its nature and characteristics. The only common denominator is precipitationdeficit, which in itself may be of different magnitude per location. Overall, the seasonal rainfall benchmark of 500 mm as stipulated by Haile et al. 48 for a rough drought indicator in a season helps to present an outlook of the drought. In 2015, the entire study area fell below the 500 mm mark, while in 2017, only one tabia/kabele (village) did. Hence the 2015 drought was more widespread than that of 2017, with the majority of farmers in the lowland reporting entire crop losses. The predictors of yield and yield losses indicate that droughts coinciding with planting and maturing stages of crops are critical and inimical to crop production and agricultural sustainability and food security. Table 9. Predictors of YpHa, SCY and YLR in the 2017 drought year in the study area. ***Regression predictor is significant at the above 0.001 level. **Regression predictor is significant at the 0.01 level. *Regression predictor is significant at the 0.05 level. www.nature.com/scientificreports/ Specifically, the stability (time) dimension of food security is of concern in this study. According to the definition of Food and Agriculture Organization of the United Nations 49 the stability dimension also referred to as the time dimension of food security by Peng and Berry 50 refers the ability of a person, community or nation to withstand food system shocks arising from either natural or man-made disaster sources. The widespread crop losses witnessed in the lowland communities due to the 2015 drought indicate the susceptibility of their food systems to drought events. This condition requires the adoption of preventive adaptation measures and improved preparedness to forestall future occurrences on the impact. It is necessary to strengthen the stability dimension of food security of the vulnerable communities because of its potential effect on other levels and dimensions of food security 50 .
In both years (2015 and 2017), and for the entire study area, both annual and the kiremt seasonal rainfall were related to the crop yield reported by farmers. Crop yield is largely a function of adequate soil moisture at all stages of crop growth and maturity. Hence, the shortage of rainfall (occurrence of drought) brought about reduced crop yield in the area. These assertions are corroborated by the findings of Gidey et al. 17 who reported that agricultural drought responds positively to seasonal rainfall. Also, part of the findings of previous studies [51][52][53] indicate that rainfall significantly determines the occurrence of drought, hence manifesting high influence on the growth and development of vegetation.
Different timescales of SPI, especially the 3-and 12-month scales calculated in this study were recurring predictors of crop yield in the area. Although, the 6-month scale is excluded in this study, we can infer that it is likely to be a predictor of drought. In relating meteorological drought index for drought monitoring, several studies have found a significant relationship between SPI and crop yield [54][55][56] . Also, the vegetation anomaly (Dev-NDVI) can detect agricultural drought and predict crop yield in the area. These rainfall and vegetation based indices used in this study considered vis-à-vis specific crop timelines are helpful to clarify the periods drought were more critical. For example, the SPI-3 for July 2015 predicted sorghum yield for the lowland portion of the study area, because sorghum crops are at their early stages of growth at this particular time. In the mid and highland areas, the Dev-NDVI for October 2015 predicted yield and yield loss as the late maturity period for barley includes October. Hence, the satellites detected low vegetation vigor, indicating drought at the time. For the Dev-NDVI to detect agricultural drought in October, it also implies that preceding moisture deficit has occurred.
The reduced impact of drought on crop yield in 2017, especially in the lowland area is attributable to the increased rainfall received. Most of the tabias (villages) received rainfall higher than 500 mm, which supports healthy and sustainable agriculture. Additionally, the report of Madajewicz et al. 57 contains Oxfam's contribution to the minimization of drought risks by the construction of canals covering 862 hectares of land to supplement rain-fed agriculture in the occurrence of a drought event, pointing to the increased popularity of spate irrigation systems. This may likely have been utilised in these areas resulting in reversing the report of 100% crop losses witnessed in 2015, and obtaining above-normal crop yield in 2017, despite droughts. Steenbergen et al. 58 describe the dual benefit of spate irrigation as being a means of both soil and water conservation. Hence, while water channeled from the highlands to the midlands and lowlands are stored for use as required, the alluvial deposit that accompanies the flood are rich in nutrients to improve crop yields and reduce costs of purchasing fertilizers 58 .

Conclusion
This study critically assessed the recent drought years (2015 and 2017) in northern Ethiopia, with a focus on a small area containing high, mid and lowland agroecology. Results obtained show droughts as having spatiotemporal variations and impacts, with its primary and common reflection being reduced vegetation amount. This translated to crop failures, the reduced income of smallholders, food shortages and food insecurity. SPI, NDVI and Dev-NDVI were found quite useful as indicators of drought and considered able to be applied in the preparation of early warning for droughts when they fall below a certain established threshold for the area in each cropping season, especially at critical crop growth and maturity stages. Findings of this study can provide information for further discussions in improving drought monitoring at lower locational scales. The stability (time) dimension of food security was affected by the 2015 drought in the lowland areas as most farmers report entire sorghum losses. Irrigation channels to flow from the highland areas to the lowland areas (spate irrigation) could be popularized as an adaptation measure to reduce vulnerability to droughts. This will improve food security/stability in the face of climatic events such as droughts by providing supplementary water supply to meet the crops' water requirements especially during the growing and maturity periods. With the provision of supplementary water source(s) during rainfall shortages in cropping seasons, the impacts of future droughts on food security are reduced. There could be more factors responsible for crop losses not captured in this study, hence providing an aspect for future further research.

Limitations and further research
The findings of this study would have been more generalizable if the number of farmers that constituted the sample size was more in the number. The complexities and systematic measurement errors in self-reported crop production data from smallholders in Africa identified by past studies [59][60][61] are hereby acknowledged as limiting the findings of this study. The exclusion of temperature as a variable in this study is an obvious limitation, which would have provided improved outcome for drought assessment, especially with the use of SPEI, instead of SPI. Also, the estimation of a 6-month time-step for SPI (SPI-6) was omitted from the study, which could further strengthen findings. However, future studies can address these deficiencies, in addition to conducting a more detailed analysis of specific crops' growth stages with climatic variables, and the impact of the use of spate irrigation systems, for more robust outcomes. Spate irrigation was not considered a priori as a factor that could affect the impacts of droughts in the area, hence further inquiry was not considered in this regard.

Data availability
The datasets used for this study are available and can be provided by the corresponding author on request.